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Current theoretical models of Saturn's ionosphere are similar to those of 

Jupiter's because of the gross similarity of their upper atmospheres. In the 

case of Jupiter, the theoretical models can be fit reasonably well to 

[121 

ionospheric electron density profiles * obtained from the Pioneer and 

[3 4 5 ] 

Voyager Radio Science occultation experiments. ’ In contrast, the 
theoretical models of Saturn's ionosphere are inconsistent with both the 
ionospheric electron density profiles 1 *' J obtained from the Pioneer and 
Voyager occultation observations^’ 10 ^ and the large diurnal variation of 
maximum ionospheric electron density deduced from studies of Saturn lightning 

[11 19 1 

discharges. * We propose a radically different model of Saturn's 

ionosphere in which water plays a major role as a minor constituent present by 

virtue of downward diffusion from an external source. Our model Saturn 

ionosphere is a classical "F2" type layer resulting from the photodissociative 

production of H + from and rapid chemical loss due to a series of charge 

7 -2 

exchange reactions with. water. A planetwide influx of ^xIO molecules cm 

s' 1 of water from the rings is consistent with the observed ionospheric 

electron densities. An enhanced influx of water occurs at latitudes (-38°, 

[ 12 ] 

+44°) magnetically connected to the inner edge of Saturn's B ring which 

results from an electromagnetic erosion process^ ^ contributing substantially 

to the (local) upper atmosphere water content. Present-day influx at these 

9 -2 -1 

latitudes may be as large as ^2x10 molecules cm s . 
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The essential features of previous theoretical models of Saturn's 

ionosphere are summarized below. The major topside ion, H + , is produced via 

photoioniziation of atomic H or photodissociative ionization of Hg. Both H 

and H + are extremely long-lived above the homopause ( 100 km altitude; all 

altitudes referenced to the 1 bar pressure level) and are removed by downward 

diffusion to the vicinity of the homopause where appreciable concentrations of 

heavier gases (most importantly methane) are encountered. H + ion loss occurs 

essentially at the homopause via charge exchange and electron-ion 

5 -3 

recombinations, resulting in a peak ion concentration of ^2x10 cm at </, 1200 

km altitude. This very simplified picture emerges from a consideration of 

many possible chemical reactions and detailed calculations taking into account 

the production, diffusion and chemical losses of many atmospheric 
[7] 

constituents. 

Figure 1 is a representative summary of the dawn and dusk electron density 
(as a function of altitude) profiles obtained from analyses of the Pioneer and 
Voyager radio occultation data. In Figure 1, mid-latitude Voyager observations 
are compared with a typical model ionospheric electron density profile. 

(Figure 2 shows all of the radio science observations.) The measured electron 
densities are all approximately an order of magnitude less than that calculated 

on the basis of current models. Maximum observed ionospheric densities are 

4_o 5-3 

^1x10 cm , compared with a calculated maximum density of 2x10 cm in all 

previous models. The peak electron density is observed at between 2000 and 

3000 km altitude instead of the calculated 1200 km altitude. The discrepancy 

in the altitude of the peak is particularly significant in light of the thermal 

structure of the atmosphere. Above >''1600 km, the atmosphere is isothermal with 

a temperature of ^800 °K and an atmospheric scale height of 
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^400 km. Below that the temperature decreases linearly to v‘125°K at 1000 km 

[ 7 ] 

altitude where the atmospheric scale height is s60 km. The difference in 
the observed and modeled altitudes of peak ion density is thus many 
atmospheric scale heights. The large (^2 order of magnitude) diurnal 
variation of maximum ionospheric electron density (Figure 3) deduced from 
analyses of Saturn lightning discharges^ 1 ^ is even more difficult to 
understand on the basis of current models of Saturn's ionosphere. Without an 
additional loss mechanism for rapid removal of the long-lived H + ions, no 
diurnal variation is expected since the time constant for loss of H + above the 
homopause is ■/'400 hrs, many planetary rotations. This is due to the slow 
electron recombination reaction H + + e -► H + hv @ 2x1 o” cm s and the 
large diffusion time constant below the 2500 km altitude of the ionospheric 
peak. 


The observational results thus require an additional and effective loss 

mechanism for H + ions, since ion production via photodissociative ionization 

of and photoionization of atomic H is well understood. Several possible 

loss processes have been quantitatively evaluated but none have proven 

adequate. These include losses due to an increased methane abundance which 

requires an unreasonably high eddy diffusion coefficient^^, removal of H + by 

vibrationally excited and/or the presence of large vertical 

[7 8] 

drifts. * The possibility of an interaction between the ionosphere and 
ring system was recognized as early as 1975^^ and a suggestion involving H + 
ion loss due to the presence of OH radicals from the rings^^’ 1 ^ was made to 
explain the low ionospheric densities inferred from the Pioneer observations. 
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We propose a simple model ionosphere dominated by the removal of H + ions 
via charge exchange with water and subsequent electron recombination. • Table 1 
summarizes the important chemical reactions and reaction rates; the reactions 
involving water have been largely drawn from work relating to comet 
models. ^ The important reactions are H + ion production via dissociative 
photoionization of H 2 (A3) and loss via the charge exchange with H 2 o (B1). 

The H 2 0 + ion is very short lived, reacting with H 2 (B2) to produce H^O*, which 
is then readily removed by electron recombination (B3-B5). As in previous 
models, the H* ions produced by photo ionization of H 2 (A2) are rapidly removed 
via reaction with H 2 (B6) and charge exchange (B7) and are not a major topside 
ion. The direct photoionization of atomic hydrogen (A4) does not 
significantly contribute to the H + ion density due to the limited availability 
of atomic H. 

The OH produced in B4 and B5 has a short lifetime (t < Is) at the 
altitude of the peak ion density due to a reaction (B9) with H 2 , producing H 2 0 
and H. Thus water acts as a catalyst in the rapid removal of H + ; each H + ion 
removed returns, on average, 3 or 4 H atoms with no associated chemical loss 
of H 2 0. Photodissociation of HgO (A5) yields OH which in the presence of 
appreciable H 2 (B9) rapidly reforms H 2 0 (as would OH from any other source). 
There is no net loss of water due to charge exchange or photoionization. Thus 
the density of water above the homopause is given by^ 20 ^ 

N = N„ EXP(-9Z/H) + FH n /8D n [EXP(-Z/H)-EXP(-9Z/H)] 
w u u u 

where N Q , H Q , and D Q are the H 2 0 density, the atmospheric scale height, and 
molecular diffusion coefficient at a reference (Z=0) altitude, F is the 
downward flux of H 2 0 and H is the atmospheric scale height. Assuming that the 
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water is readily removed by freezing out just above the homopause or by mixing 
and chemical loss at the homopause, the density of water throughout the 
(isothermal) region of interest above 1600 km altitude is just N w = 
(FH/8 Dq)EXP(-Z/H) . Since the diffusion coefficient D=Dq/N, where N is the 
atmospheric density, the downward flux of water results in a water density 
proportional to the density of the H 2 atmosphere through which it flows under 
the condition of constant flux, i.e., N w /N = cnst = FH/8 NqD q . 

Under these conditions a classical "Bradbury Layer", analogous to the 

earth's F2 layer can form at an altitude where the chemical and diffusive 

[20 21 ] 

time constants are comparable. ’ For such a layer, the ionospheric 

peak forms above the level of maximum ion production due to the upward 

diffusion of ions. Below that altitude the increasing chemical loss (charge 

transfer with water) effectively removes the major ion, so that most of the 

ionization anticipated by previous models is lost. The loss coefficient is 6 

= ktHgO] where k is the rate coefficient of reaction B1 in Table 1. The 

required water concentration at the altitude of the ionospheric peak can be 

estimated by N = q(h )/S(h ) where N is the maximum electron density at the 
m m m m 

altitudes (h m ) of the ionospheric peak, and the ion production function q and 

loss g are evaluated at the altitude of the peak. Taking h m *2500 km, q s 1 
-3-1 5-3 

cm s and >H0 cm appropriate to the equatorial ionosphere at noon 

-3 

yields a water concentration of ^1250 cm at 2500 km altitude. This requires 
7 -2 -1 

a 4x10 cm s downward flux of water molecules. The maximum densities 

observed at the dawn and dusk terminators where the EUV production rate is 

J’0. 1 of the maximum are -'’I order of magnitude lower than the noon density, as 

2 

expected. Equating the chemical and diffusive time constants (1/S </* H /D) 
yields an expected peak altitude between 2,000 and 2500 km. 
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The variation of electron density with altitude can also be understood 

using simple approximations to "Bradbury layer" behavior that have been widely 

[20 21 ] 

used in the study of the earth’s F2 layer. ' The ion distribution well 
above the peak is controlled by ambipolar diffusion, so N^expC-Z/Hp) where Hp 
is the plasma scale height. This distribution is indicated in Figure 2 
assuming 800 °K electron and ion temperatures, for which Hp = 4H. The H + ion 
density distribution below the peak is controlled by chemical loss, with N(h) 
^q/S(h). The ion production, q, varies slowly, so to first approximation we 
expect N(h) below the peak to increase with height as exp(h/H), as is 
illustrated in Figure 2. This reflects the increased loss at greater depths 
due to greater water concentration. The peak has a width of a few scale 
heights and occurs at an altitude (2500 km) where the diffusion time constant 
is a few hours, as does the earth's F2 layer. However, the rapid nighttime 
decay of the equatorial ionospheric density evidenced in Figure 3 requires a 
more rapid loss, perhaps a downward drift related to the 
ionosphere/protonosphere exchange described below. 

Figure 3 shows an ionospheric density minimum just after dusk and a 

steady increase in near-equatorial ionospheric density throughout the night. 

We attribute this phenomenon to a diurnal ionosphere - protonosphere exchange 

[ 20 ] 

similar to that proposed for the earth's nighttime F2 layer. During 

daylight hours, photoproduction in Saturn's ionosphere pumps protons outward 
along field lines to great heights. At night, a relaxation back towards the 
ionosphere occurs to replenish the nighttime ionospheric density. If we 
assume a uniform density (n) column of protons at rest at dusk, under constant 
downward acceleration (a) the proton influx # = nat increases linearly with 

time. Supposing the nighttime maximum ionospheric density occurs at high 
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altitudes where loss is negligible, the increase in the maximum density with 

time is related to the influx by N (t) = /H dt = (na/2H)t*. This increase 

m p 

is indicated on figure 2 by the dashed line, which fits the nighttime obser- 

-5 -3 2 

vations well with na/2H ^2x10 cm s , corresponding to a proton influx 

7 -2 -1 

increasing to a maximum of ^2x10 cm s just before dawn. The equivalent 
ion production rate at dawn due to this exchange (v>.2 cm ^ s S is comparable 
to or greater than the EUV production. Thus, we envision a very dynamic 
ionosphere, characterized by ion production and transport into a vast 
protonosphere reservoir by day coupled with relaxation and return by night. 

The increased ion losses at latitudes magnetically connected to the inner 
[ 12 ] 

edge of Saturn's B-ring can be explained by an increased local water 

9 - 2—1 

influx of as much as 2x10' molecules cm s . The erosional process 

described by Northrop and Hill^^ provides a mechanism whereby sub-micron 

charged particles (collection of water molecules), unstable in Saturn's ring 

plane, precipitate along magnetic field lines into the ionosphere. If these 

particles are deposited directly at the level of the ionospheric peak (<^2500 

9 -2 -1 

km) the required flux is less than 2x10 cm s . At this maximum flux, the 

—14 —2 -1 

ring mass loss is ^6x10 g cm s . The large dayside extinction was 

[ 12 ] 

observed at latitudes of -37 to -39,° the maximum southerly extent of the 

Voyager 1 spacecraft. Averaged over 4° latitude extent of the presumed 

source region in the rings (1.52 - 1.62 Rg), the total ring mass loss is 
6 — 1 12 —1 

^1.2x10 g s~ or 1.6x10 g yr , comparable to the estimated total mass 

[ 22 ] -2 

influx from interplanetary meteoroids. At this rate, 60 g cm of 

material can be removed from the source region in ^30 10^ years, suggesting 
that this process may be responsible for the large decrease in optical depth 
of the rings inward of 1.62 R,,. 
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The radio science observations of Saturn's ionosphere and the inferences 
drawn from studies of Saturn lightning can be largely understood on the basis 

of a new model of Saturn's ionosphere in which water plays a major role. The 

7 -2 -1 

planetwide influx of water (4x10 molecules cm s ) required is 

approximately an order of magnitude greater than that expected from 

[23] [22] 

photosputter ing of the rings. Morfill et al. estimate that 

micrometeoride impacts on Saturn's icy rings produce between .75 to 6x10^ 

g/sec of vapor, <^10$ of which is deposited on Saturn's atmosphere. This 

7 -2 -1 

results in an estimated influx of between 0.6 and 4.5x10 molecules cm s 

[ 22 ] 

In spite of the uncertainties in this estimate, micrometeoride bombardment 

of Saturn's rings is an attractive source mechanism because it provides an 

explanation of the vast differences between Jupiter's and Saturn's 

ionospheres. An atmospheric influx of water has been previously suggested to 

[24 ] 

explain features of ultraviolet spectra of Saturn and Voyager infrared 
[251 

spectra of Titan. The increased influx of water at latitudes magnetically 

connected to the inner edge of the B-ring is not inconsistent with expected 
erosion rates. Finally, we infer a diurnal ionosphere - protonosphere 
exchange in which protons pumped into the protonosphere by day return by night. 
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TABLE 1 


CHEMICAL REACTIONS IN SATURN'S IONOSPHERE : Hg, HgO, H 

ION PRODUCTION RATE* REF. 


A1 

Hg + hv ♦ H + H 

in 

• 

X 

o 

i 

o 
c a 

L 

C7 3 

A2 

H* + e 

5.4 

X 

10' 10 

[7] 

A3 

H + H + + e 

9.5 

X 

10" 11 

[7] 

A4 

H + hv ♦ H + + e 

7.3 

X 

10- 10 

[7] 

A5 

H 2 0 + hv * OH + H 

1.0 

X 

10- 7 

[193 


0H + + H + e 

5.5 

X 

10- 10 

[193 


H 2 + 0( l D) 

1.4 

X 

10" 8 

[193 


H 2 0 + + e 

3.3 

X 

10“ 9 

[193 


H 2 + 0 + + e 

5.8 

X 

10" 11 

[193 


H + + OH + e 

1.3 

X 

10- 10 

[193 


CHARGE EXCHANGE/RECOMBINATION 


B1 

H_0 • 

+ H+ + H_0 + + H 

8.2 

X 

10" 9 cm 8 s" 1 

[183 

B2 

t— 

H 2 0 + 

C. 

+ H 2 + H 3 0 + + H 

6.1 

X 

10- 10 

T 

[183 

B3 

B4 

H 3 0 + 

+ e ♦ H 2 0 + H 
H 2 + OH 

4.3 

4.3 

X 

X 

o o 
i -tj -! 

[183 

[183 

B5 


H + H + OH 

4.3 

X 

10" 7 

[183 

B6 

h 2 + 

H 2 * H 3 + + H 

2.0 

X 

10 -9 

[73 

B7 

H 3 + 

e + H 2 + H 

1.0 

X 

10 -7 

C73 

B8 

H* + 

e -► H + v 

2.0 

X 

10" 12 

[7 3 

B9 

OH + 

h 2 * h 2 o + h 

3.6 

X 

10’ 11 

[183 


* Photoionization rates normalized to 10 A.U. by 1/R 2 scaling of reference 
values 
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Figure Captions 


Figure 1 Model calculations of mid-latitude ionospheric electron density as a 
function of altitude (EUV) and Voyager observations of Saturn’s 

[ 7 ] 

ionsphere at mid-latitudes, adapted from Atreya et al. 

Figure 2 Ionospheric electron density profiles deduced from Voyager 1 and 2 

[7] 

and Pioneer radio science data, adapted from Atreya et al. 

Altitude is referenced to the 1 bar pressure level. 

Figure 3 Diurnal variation of maximum ionospheric electron density inferred 

from Voyager observations of lightning-generated radio emissions in 

[ 12 ] 

Saturn's atmosphere, adapted from Kaiser et al. The dashed line 

indicates a model fit to the nighttime increase of maximum 
ionospheric electron density. 
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